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Abstract. Nuclei exhibit quantum phase transitions (earlier called ground state phase 
transitions) between different shapes as the number of nucleons is modified, resulting in 
£h ' changes in the ground and low lying nuclear states. Special solutions of the Bohr Hamiltonian 

appropriate for the critical point of such shape phase transitions, as well as other special 
solutions applicable to relevant nuclear regions are described. 
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£ 1 INTRODUCTION 

in ' 

Atomic nuclei are known to exhibit changes of their energy levels and elec- 
tromagnetic transition rates among them when the number of protons and/or 
neutrons is modified, resulting in shape phase transitions from one kind of 
collective behaviour to another. These transitions are not phase transitions 
of the usual thermodynamic type. They are quantum phase transitions pQ 
(initially called ground state phase transitions [2j), occurring in Hamiltonians 
of the type H = c{H\ + gH2), where c is a scale factor, g is the control param- 
eter, and H\, Hi describe two different phases of the system. The expectation 
value of a suitably chosen operator, characterizing the state of the system, is 
used as the order parameter. 
^ ' In the framework of the Interacting Boson Model [SJ, which describes 

nuclear structure of even-even nuclei within the U(6) symmetry, possessing 
the U(5), SU(3), and 0(6) limiting dynamical symmetries, appropriate for 
vibrational, axially deformed, and 7-unstable nuclei respectively, shape phase 
transitions have been studied 25 years ago [2] using the classical limit of the 
model [UEJE], pointing out that there is (in the usual Ehrenfest classification) 
a second order shape phase transition between U(5) and 0(6), a first order 
shape phase transition between U(5) and SU(3), and no shape phase transition 
between 0(6) and SU(3). It is instructive to place [1] these shape phase 
transitions on the symmetry triangle of the IBM [7J, at the three corners of 
which the three limiting symmetries of the IBM appear. 
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More recently it has been realized [U [9] that the properties of nuclei 
lying at the critical point of a shape phase transition can be described by 
appropriate special solutions of the Bohr Hamiltonian [10] . labelled as critical 
point symmetries. The E(5) critical point symmetry [8] has been found to 
correspond to the second order critical point between U(5) and 0(6), while 
the X(5) critical point symmetry [9] has been found to correspond to the first 
order transition between U(5) and SU(3). 

The introduction of the critical point symmetries E(5) [8] and X(5) [9] 
has triggered much work on special solutions of the Bohr Hamiltonian, cor- 
responding to different physical situations. Several of these solutions will be 
mentioned here, together with experimental examples appropriate for each 
case. 



2 E(5) AND RELATED SOLUTIONS 



The original Bohr Hamiltonian [10] is 
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where (3 and 7 are the usual collective coordinates describing the shape of the 
nuclear surface, Qk (fc = 1, 2, 3) are the components of angular momentum, 
and B is the mass parameter. 

It has been known for a long time [TT] that exact separation of vari- 
ables occurs in the corresponding Schrodinger equation if potentials of the 
form V((3, 7) = U(/3) are used, i.e. if the potential is independent of the vari- 
able 7 , thus corresponding to 7-soft nuclei. Then wavefunctions of the form 
\I/(/3, 7 , 9i) = /(/3)3>(7, ^) are used, where 9i (i = 1,2,3) are the Euler angles 
describing the orientation of the nucleus in space. 

In the equation involving the angles, the eigenvalues of the second order 
Casimir operator of SO(5) occur, having the form A = r(r + 3), where r = 0, 
1, 2, . . .is the quantum number characterizing the irreducible representations 
(irreps) of SO (5), called the "seniority" [12J. This equation has been solved 
by Bes p3]. 

The values of angular momentum L contained in each irrep of SO (5) (i.e. 
for each value of r) are given by the algorithm [3j r = 3z/a + A, where = 0, 
1, ...is the missing quantum number in the reduction SO(5) D SO(3), and 
L = A, A + 1, . . . , 2A — 2, 2A (with 2A — 1 missing). The values of L allowed for 
each (t,va) have been tabulated in [3| [HI [T5]. 
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The "radial" equation can be simplified by introducing [8j reduced ener- 
gies e = Tpr-E and reduced potentials u = jpU. The form of the solution of 
the radial equation depends on the choice made for U{f5). 

2.1 E(5) 

In the case of E(5) [15], a 5-dimensional (5-D) infinite well [u((3) = if 
(3 < 0w, u{(3) = oo for > flw] is used, since the potential is expected to 
be flat at the point of a second order shape phase transition. Then the (3- 
equation becomes a Bessel equation of order v = r + 3/2, with eigenfunctions 
proportional to the Bessel functions J T +3/2{ z ) (with z = (3k, k = -y/e), while 
the spectrum is determined by the zeros of the Bessel functions 

where x^ T is the £-th zero of the Bessel function J T +3/2(z)- The spectrum 
is parameter free, up to an overall scale factor, which is fixed by normalizing 
the energies to the excitation energy of the first excited 2 + state, 2f. The 
R A/2 = E(4f)/E(2f) ratio turns out to be 2.199 . The same holds for the 
B(E2) values, which are normalized to the B{E2) connecting the two lowest 
states, B(E2;2f — > 0^). The symmetry present in this case is [15] E(5) D 
SO(5) D SO(3) D SO(2). 

2.2 Other solutions 

For u(P) = P 2 /2 one obtains the original solution of Bohr \W\ 116]. which cor- 
responds to a 5-D harmonic oscillator characterized by the symmetry U(5) D 
SO(5) D SO(3) DSO(2) [T7j, the eigenfunctions being proportional to Laguerre 
polynomials [18j, and the spectrum having the simple form En = N + 5/2, 
with N = 2v + r, and v = 0, 1, 2, 3, . . ., which has -R4/2 = 2. The spectra of 
the u{j3) = P 2 /2 potential and of the E(5) model become directly comparable 
by establishing the formal correspondence v = £ — 1. 

The Davidson potential u((3) = (3 2 + (where (5q is the position of the 
minimum of the potential) [191 [2T?1 [2"T] also leads to eigenfunctions which are 
Laguerre polynomials, the energy eigenvalues being |2CH I2T] (in Tluj = 1 units) 



E„ T = 2n + 1 + 



r + f) 2 + A* 



1/2 

(3) 



For (3q = the above mentioned original solution of Bohr [U(5)] is obtained, 
while for (3$ — > 00 the 0(6) limit of the IBM is obtained |20j . Therefore the 
Davidson potential provides a one-parameter bridge between U(5) and 0(6). 
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One can exploit this fact, by introducing a variational procedure \22\ 123]. in 
which the rates of change of the Rl = E(Lf)/E(2f) energy ratios of the 
ground state band with respect to the parameter /3q are maximized for each L 
separately. The results lead to an energy spectrum very close to that of E(5) 
[22j . The method has also been applied to other bands, as well as to B(E2) 
transition rates 1 23 1 . 



The sequence of potentials U2n{P) = 2~ (with n being an integer) leads 
for n = 1 to the Bohr case, while for n — > oo leads to the infinite well of 
E(5) [24] . Therefore this sequence of potentials provides a "bridge" between 
the U(5) symmetry and the E(5) model, using their common SO(5)dSO(3) 
chain of subalgebras for the classification of the spectra. Solutions for n / 1 
have been obtained numerically [25} [26| I14j . The solutions for n = 2, 3, 4, 
labelled as E(5)-/3 4 , E(5)-/3 6 , and E(5)-/3 8 , lead to R 4/2 = 2.093, 2.135, and 
2.157 respectively. Complete level schemes have been given in Ref. |14j . 

A bridge complementary to the one just mentioned, i.e. a bridge spanning 
the region between E(5) and the 7-soft rotor 0(5), has been obtained by using 
an infinite well potential with boundaries (3m > (3 m > |27j . The model, 
called 0(5)-CBS, since it is a 7-soft analog of the confined /3-soft (CBS) rotor 
model [281 EH]> contains one free parameter, rg = (3 m / '(3m, the value rg = 
corresponding to the E(5) model, and 773 — ► 1 giving the 7-soft rotor [0(5)] 
limit. 

Other solutions obtained in this framework are listed below. 

a) A version of E(5) using a well of finite depth, instead of an infinite 
one, has been developed [30J. 

b) The sextic oscillator, which is a quasi-exactly soluble [31, 32j potential, 
has also been used as a 7-independent potential [33J 




c) Coulomb-like and Kratzer-like potentials have been used in Ref. [34j . 

d) A linear potential has been considered in Ref. [35], where a review of 
potentials used in this framework is given. 

e) A hybrid model employing a harmonic oscillator for L < 2 and an 
infinite square well potential for L > 4 has been developed [36]. 

2.3 Experimental manifestations of E(5) 

The first nucleus to be identified as exhibiting E(5) behaviour was 134 Ba [37J, 
while 102 Pd [38J also seems to provide a very good candidate. Further studies 




on 134 Ba |39] and 102 Pd [40j . in which no backbending occurs in the ground 
state band, which remains in excellent agreement with the parameter-free E(5) 
predictions up to high angular momenta, reinforced this conclusion. 104 Ru 
[41] . 108 Pd 02], 114 Cd [33], and 130 Xe [SJ have also been suggested as possible 
candidates. A systematic search [45 \ 146] on available data on energy levels and 
B(E2) transition rates suggested 102 Pd, 106 ' 108 Cd, 124 Te, 128 Xe, and 134 Ba as 




>2n 
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possible candidates, singling out Xe as the best one, in addition to Ba. 
This is in agreement with a recent report [17] on measurements of El and Ml 
strengths of 124_136 Xe carried out at Stuttgart, which provides evidence for a 
shape phase transition around A ~ 130. 128 Xe has been measured (November 
2006) in Jyvaskyla [H]. Recently, 58 Cr has also been suggested as a candidate 

02]. 

The assumption of a flat /3-potential in the E(5) symmetry has been 
tested by constructing potential energy surfaces (PESs) for nuclei close to the 
E(5) symmetry, through the use of relativistic mean field theory [SO]. It has 
been found that the relevant PESs come out quite flat, corroborating the E(5) 
assumption. 

2.4 Odd nuclei: E(5/4) and E(5/12) 

The models discussed so far are appropriate for even-even nuclei. Odd nuclei 
can be treated by coupling E(5), describing the even-even part of an odd 
nucleus, to the odd nucleon by the five-dimensional spin-orbit interaction 
[5H OS]- If the odd nucleon is in a j = 3/2 level, the E(5/4) model [5H Qj)] 
occurs, while if the odd nucleon lives in a system of levels with j = 1/2, 
3/2, 5/2, the E(5/12) model [52] is obtained. Shape phase transitions from 
spherical to 7-unstable shapes in odd nuclei have also been considered [53l [M] 
in the framework of the Interacting Boson Fermion Model |55[ [56] for the case 
of an odd nucleon in a j = 3/2 level. Shape phase transitions in odd nuclei 
have also been considered [57] for the case of an odd nucleon in a system of 
levels with j = 1/2, 3/2, 5/2 in the framework of the U(5/12) supersymmetry 
[551 156] , giving good results in the Os-Hg region. 

A first effort to locate nuclei exhibiting the E(5/4) symmetry has been 
carried out for 135 Ba [58], with mixed results. The Ir-Au region might be a 
more appropriate one, since the U(6/4) supersymmetry has been found there 
[SUES]. 63 Cu, despite its small size, could also be a good candidate for E(5/4), 
as discussed in Ref. [5T] . 

3 X(5) AND RELATED SOLUTIONS 

In the case of X(5) one tries to solve the Bohr Hamiltonian of Eq. ([T]) for 
potentials of the form u(f3, 7) = u((3) + 14(7), seeking solutions of the relevant 
Schrodinger equation having the form ^(/3, 7, 0j) = 4>^(P,j)T>^ IK (di), where 
9i (i = 1, 2, 3) are the Euler angles, T>{9j) denote Wigner functions of them, L 
are the eigenvalues of angular momentum, while M and K are the eigenvalues 
of the projections of angular momentum on the laboratory-fixed z-axis and the 
body-fixed z'-axis respectively. One is interested in cases near axial symmetry, 
i.e. close to 7 = 0. Thus one uses a harmonic oscillator potential u{"f) = 
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(3c) 2 7 2 /2. Near 7 = the last term in the Bohr Hamiltonian can be rewritten 
in the form [9] 



E 2 / Qk 27r , x *> + Ql + Ql) + Q§ (7V - 1) ■ (4) 

fc=7^,3 sm 2 (7 " f *:) 3 Vsm 7 37 

Using this result in the Schrodinger equation corresponding to the Hamiltonian 
of Eq. ([1]), introducing reduced energies e = 2BE/H 2 and reduced potentials 
u = 2BV/h 2 as in the E(5) case, and taking into account that the reduced 
potential is of the form u{(3,^f) = u{0) + u(j), the Schrodinger equation can 
be separated into two equations [D]. 

In the equation containing the 7- variable, (3 2 denominators remain, which 
are replaced by their average values over the (3 wavef unctions, (f3 2 )- Taking 
into account the simplifications imposed by 7 being close to zero, the rel- 
evant equation takes the form corresponding to a two-dimensional harmonic 
oscillator in 7, having wavefunctions proportional to Laguerre polynomials [9]. 

The form of the solution of the radial equation depends on the choice 
made for u(f3). 

3.1 X(5) 

In X(5) [9] a 5-D infinite well potential is used, as in E(5). The relevant 
equation is again a Bessel equation, but with order 

»=(^ + lf. 

The solutions still have the form of Eq. ([2]), with (£, r) replaced by (s,L), 
where s is the order of the relevant root of the Bessel function J v {k s ,Ll3)- The 
relevant exactly soluble model is labelled as X(5) (which is not meant as a 
group label, although there is relation to projective representations of E(5), 
the Euclidean group in 5 dimensions [9j). The total energy has the form 

E(s, L, n 7 , K, M) = E + B{x S)L f + An^ + CK 2 , (6) 

where n 7 is the quantum number of the two-dimensional oscillator occurring 
in the 7-equation, while Eq, A, B, C are free parameters. From this equation 
it is clear that the spectra of the ground state and f3 bands only depend on 
an arbitrary scale, fixed by normalizing them to the energy of the 2f state, 
while the bandheads of the 7 bands are parameter dependent. (The spacings 
within the 7 bands are however fixed [59J.) The R4/2 ratio is 2.904 . 
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3.2 Other solutions 



For u((3) = f3 2 /2 one obtains an exactly soluble model which has been called 
X(5)-/3 2 [60J, the eigenfunctions being proportional to Laguerre polynomials 
and the spectrum having the form 



with R4/2 = 2.646. The spectra of the u(f3) = (3 2 /2 potential and of the X(5) 
model become directly comparable by establishing the formal correspondence 
n = s — 1, where n is the usual oscillator quantum number. 

The Davidson potential u(/3) = 1 + % (where (5q is the position of the 
minimum of the potential) [X9^ |2"0"1 I2T] also leads to eigenfunctions which are 
Laguerre polynomials, the energy eigenvalues being [22J 123] (in hu = 1 units) 



For 0o = the above mentioned X(5)-/3 2 solution is obtained, while for /3o — > 
00 the rigid rotor limit is obtained. Therefore the Davidson potential provides 
a one-parameter bridge between X(5)-/3 2 and the rigid rotor. One can exploit 
this fact, by introducing a variational procedure [22], [23] , in which the rates of 
change of the Rl = E(Lf)/E(2f) energy ratios of the ground state band with 
respect to the parameter (3q are maximized for each L separately. The results 
lead to an energy spectrum very close to that of X(5) [22]. The method has 
also been applied to other bands, as well as to B{E2) transition rates [23] . 

The sequence of potentials U2 n (0) = (with n being an integer) leads 
for n = 1 to the X(5)-/3 2 case, while for n — > 00 leads to the infinite well of 
X(5) [21]. Therefore this sequence of potentials provides a "bridge" between 
the X(5)-/3 2 solution and the X(5) model, in the region lying between U(5) 
and X(5). Solutions for n 7^ 1 have been obtained numerically [60]. The 
solutions for n = 2, 3, 4, labelled as X(5)-/3 4 , X(5)-/? 6 , and X(5)-/? 8 , lead to 
-R4/2 = 2.769, 2.824, and 2.852 respectively. Complete level schemes have been 
given in Ref. [60j . 

A bridge complementary to the one just mentioned, i.e. a bridge span- 
ning the region between X(5) and the rigid rotor, has been obtained by using 
an infinite well potential with boundaries (3m > An > [28, EH]. The model, 
called the confined /3-soft (CBS) rotor model [2E1I2S]) contains one free pa- 
rameter, rg = (3 m /(3M, the value rp = corresponding to the X(5) model, and 
rp — > 1 giving the rigid rotor limit. 

Other solutions obtained in this framework are listed below. 





(7) 



1 1/2 



E n ,L = 2n + 1 + -L(L + 1) + - + /? 4 



(8) 



7 



a) A potential with linear sloped walls has been considered in Ref. |61j . 
The sloped walls result in a slower increase of the energy levels of the /3 band 
as a function of L in this model, as compared to X(5). This feature improves 
agreement to experiment. 

b) Coulomb-like and Kratzer-like potentials have been used in Ref. |62j . 

c) The approximate separation of variables used in X(5) has been tested 
recently through exact numerical diagonalization of the Bohr Hamiltonian [63j , 
using a recently introduced [Ml ESJ [66] computationally tractable version of 
the Bohr-Mottelson collective model. 

d) Exact separation of variables can be achieved by using potentials of 
the form u(f3, 7) = u{0) +u( , y)/(3 2 [U\. This possibility has been recently 
exploited for the construction of exactly separable (ES) analogues of the X(5) 
and X(5)-/? 2 models, labelled as ES-X(5) and ES-X(5)-/3 2 respectively [67]. 

e) The use of periodic u(j) potentials in the X(5) framework has been 
recently considered in Ref. [68] . 

A review of potentials used in this framework is given in Ref. [35j . 

3.3 X(3) 

The special case in which 7 is frozen to 7 = 0, while an infinite square well 
potential is used in (5, leads to an exactly separable three-dimensional model, 
which has been called X(3) [69] . This model involves three variables, (3 and the 
two angles used in spherical coordinates, since the condition 7 = guarantees 
an axially symmetric prolate shape, for which the two angles of the spherical 
coordinates suffice for determining its orientation in space. 

Exact separation of variables is possible in this case. The equation in- 
volving the angles has the usual spherical harmonics as eigenfunctions, the 
relevant eigenvalues being L{L + 1), while the /3-equation, in which an infinite 
square well potential is used, takes the form of a Bessel equation. The radial 
solutions have the same form as in X(5), but with order 

Il(L + 1) 1 

which should be compared to Eq. ([5]). It should be noticed that in E(3), 
the Euclidean algebra in 3 dimensions, which is the semidirect sum of the T3 
algebra of translations in 3 dimensions and the SO (3) algebra of rotations in 3 
dimensions [70], the eigenvalue equation of the square of the total momentum, 
which is a second-order Casimir operator of the algebra, also leads [70 \ [Tl] to 

a similar solution, but with u = L + | = yj L(L + 1) + |. 

From the symmetry of the wave functions with respect to the plane which 
is orthogonal to the symmetry axis of the nucleus and goes through its center, 
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follows that the angular momentum L can take only even nonnegative values. 
Therefore no 7-bands appear in the model, as expected, since the 7 degree 
of freedom has been frozen. The R4/2 ratio is 2.44 . Complete level schemes 
have been given in [69j . 

3.4 Experimental manifestations of X(5) 

The first nucleus to be identified as exhibiting X(5) behaviour was 152 Sm [7T] . 
followed by 150 Nd [72]. Further work on 152 Sm [7J [7J E3 [59] and 150 Nd 
\7A \ \75 \ [76] reinforced this conclusion. The neighbouring N=90 isotones 154 Gd 
[771 ITS] and 156 Dy [78], [79] were also seen to provide good X(5) examples, the 
latter being of inferior quality. In the heavier region, 162 Yb [SO] and 166 Hf 
[81] have been considered as possible candidates. More recent experiments 
on 176 Os and 178 Os [82] indicate that the latter is a good example of X(5). 
A systematic study [83] of available experimental data on energy levels and 
B(E2) transition rates suggested 126 Ba and 130 Ce as possible good candidates, 
in addition to the N=90 isotones of Nd, Sm, Gd, and Dy. A similar study in 
lighter nuclei [M] suggested 76 Sr, 78 Sr and 80 Zr as possible candidates. 104 Mo 
has been suggested as a candidate for X(5) based on available spectra [85, 84] , 
but later studies on B{E2) values gave results close to the rigid rotor limit 
[86] . 122 Ba [57] is currently under consideration, since its ground state bands 
coincides with this of X(5). Recent measurements [88] on 128 Ce indicate that 
this nucleus is a good example of X(5). This is expected, since 128 Ce, having 
8 valence protons and 12 valence neutron holes, matches 152 Sm, possessing 12 
valence protons and 8 valence neutrons, which is a good example of X(5). 

The assumption of a flat /3-potential in the X(5) symmetry has been 
tested by constructing potential energy surfaces (PESs) for nuclei close to the 
X(5) symmetry, through the use of relativistic mean field theory [50, 89, 90, 91J. 
It has been found that the relevant PESs exhibit a bump in the middle, in 
accordance to calculations using an effective deformation, determined by 
variation after angular momentum projection and two-level mixing [92], as 
well as in Nilsson-Strutinsky-BCS calculations [93J for 152 Sm and 154 Gd. 

4 Z(5) AND RELATED MODELS 

Z(5) [94J is an analogue of X(5) appropriate for triaxial nuclei. In both cases 
potentials of the form u(/3, 7) = u((5) + u(-y) are considered. In the X(5) case 
the Hamiltonian is simplified by focusing attention near 7 = 0, which corre- 
sponds to prolate axially symmetric nuclei. In Z(5) attention is focused near 
7 = n/6, which corresponds to maximally triaxial shapes. It is known [95] 
that for 7 = 7r/6 the projection of the angular momentum on the body-fixed 
x'-axis, labelled as a, is a good quantum number, while the projection on 
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the body-fixed z'-axis, labelled as K, is not a good quantum number. One 
then seeks solutions of the relevant Schrodinger equation having the form 
tf(/3,7,0i) = ^(ftT^a^i), where 0* (i = 1, 2, 3) are the Euler angles, 
T>(9i) denote Wigner functions of them, L are the eigenvalues of angular mo- 
mentum, while M and a are the eigenvalues of the projections of angular 
momentum on the laboratory-fixed z-axis and the body-fixed x'-axis respec- 
tively, a has to be an even integer [95J. Instead of the projection a of the 
angular momentum on the x'-axis, it is customary to introduce the wobbling 
quantum number ES] n w = L — a. 

One is interested in cases near maximal triaxiality, i.e. close to 7 = 7r/6. 
Thus one uses a harmonic oscillator potential u(p/) = c (7 — /2 = c7 2 /2, 
with 7 = 7 — 7r/6. 

In the case in which the potential has a minimum around 7 = ir/6 one 
can write the last term of Eq. (pQ) in the form 

£ 2 / 9 " 2 \ ^ + 4(Q| + Q|)=4(^ + ^ + Q§)-3Qf. (10) 
fc= l i2)3 sin^ (7 - f fej 

Using this result in the Schrodinger equation corresponding to the Hamilto- 
nian of Eq. (pQ), introducing [S] reduced energies e = 2BE/H 2 and reduced 
potentials u = 2BV/h 2 , and assuming [9] that the reduced potential can be 
separated into two terms, one depending on f3 and the other depending on 7, 
i.e. u((3, 7) = u(f3) + 14(7), the Schrodinger equation can be separated into two 
equations [M] . 

In the equation containing the 7- variable, (3 2 denominators remain, which 
are replaced, in analogy with X(5), by their average values over the (3 wave- 
functions, (/3 2 ). Taking into account the simplifications imposed by 7 being 
close to vr/6, the relevant equation takes the form corresponding to a simple 
one-dimensional harmonic oscillator in 7, having wavefunctions proportional 
to Hermite polynomials [94|. 

The form of the solution of the radial equation depends on the choice 
made for u(f3). 

4.1 Z(5) 

In Z(5) a 5-D infinite well potential is used, as in X(5). The relevant equation 
is again a Bessel equation, but with order 

y/4L(L + 1) - 3a 2 + 9 _ \JL{L + 4) + 3n w (2L - n w ) + 9 
V ~ 2 ~ 2 ' 

The solutions still have the form of Eq. (J2J), with (£,r) replaced by (s,i^) = 
(s,nyy,L), where s is the order of the relevant root of the Bessel function 
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Jv{k s ,vf3)- The relevant exactly soluble model is labelled as Z(5) (which is not 
meant as a group label). 

The total energy has the form 



where is the quantum number of the one-dimensional oscillator occurring 
in the 7-equation, while Eo, A, B are free parameters. 

The wobbling quantum number n w labels a series of bands with L = 
n w ,n w + 2,n w + 4, . . . (with n w > 0) next to the ground state band (with 
iT'w = 0) [95J. The ground state band corresponds to s = 1, n w = and has 
R4/2 = 2.350 . We shall refer to the model corresponding to this solution as 
Z(5) (which is not meant as a group label), in analogy to the E(5) [BJ, and X(5) 
[9] models. Complete level schemes have been given in |94t I97j . A preliminary 
comparison to experiment has suggested 192_196 pt as possible Z(5) candidates 

(SUET]. 

4.2 Other solutions 

Solutions in the vicinity of 7 = ir/6 have also been worked out considering 
potentials of the form 7) = u{j3) + u(^)/ (5 2 [981 EH], which allow for an 
exact separation of variables. Coulomb, Kratzer, harmonic, Davidson, and 
infinite square well potentials have been considered as u(f3) in this approach 
[SHIES] ? while a displaced harmonic oscillator has been used as 14(7). A periodic 
potential 11(7) = \ij sin 2 (37) has also been considered |100j . A solution similar 
to Z(5), but with a 11(7) proportional to cos 2 (37), has also been considered 



The special case in which 7 is frozen to 7 = vr/6, while an infinite square well 
potential is used in f3, leads to an exactly separable four-dimensional model, 
which has been called Z(4) [102] . This model involves four variables, (3 and 
the three Euler angles, since 7 in this model is not treated as a variable but 
as a parameter, as in the model of Davydov and Chaban [103J . 

Exact separation of variables is possible in this case. The equation involv- 
ing the Euler angles has been solved by Meyer-ter-Vehn [95], the eigenfunc- 
tions being appropriate combinations of Wigner functions. The /5-equation, 
in which an infinite square well potential is used, takes the form of a Bessel 
equation. The radial solutions have the same form as in Z(5), but with order 



E(s, nw, L, nzf) = Eq + A(x s ^) 2 + Briy 



'7j 



(12) 



[101J. 



4.3 Z(4) 




(13) 
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which should be compared to Eq. (Ilip , where the various symbols have the 
same meaning as in Z(5). The -R4/2 ratio is 2.226 . Complete level schemes 
have been given in [102, 97J. A preliminary comparison to experiment has 
suggested 128 ~ 132 Xe as possible Z(4) candidates |102|, [9T] . 

4.4 Transition from axial to triaxial shapes 

A special solution of the Bohr Hamiltonian corresponding to a transition from 
axially deformed to triaxially deformed shapes has been given in Ref. pL04|, 
called Y(5). The proton-neutron triaxiality occurring in the SU(3)* limit |1U5| 
of IBM-2 [3] has triggered the detailed study of the phase structure of IBM-2 
[1061 11071 1108] . The main features of proton-neutron triaxiality are a low- lying 
K = 2 band and B(E2)s resembling closely these of the Davydov model [109J. 

5 CONCLUSIONS 

A great and still growing interest has been developed in the last five years 
in special solutions of the Bohr Hamiltonian, in relation to shape phase tran- 
sitions and critical point symmetries in nuclei. Extensions of these ideas in 
many directions are ongoing, including the consideration of dipole |47j and 
octupole [1101 1 1 1 1|. 1112} I113j degrees of freedom. Many developments rele- 
vant to shape phase transitions and critical point symmetries, not mentioned 
in this work (which is not a review article but rather a biased brief account 
of topics related to the authors' work), can be traced from the references in 
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